What is claimed is: 

[Claim 1] 1 . An MR system comprising: 

a magnetic resonance imaging (MRI) system having a plurality of 
gradient coils positioned about a bore of a magnet to impress a polarizing 
magnetic field and an RF transceiver system to acquire MR images; and 
a computer programmed to: 

generate phase-difference images from a first acquisition and 
a second acquisition; 

determine a non-zero background phase from the phase- 
difference images that are due to eddy currents induced by flow encoding 
gradients used to generate the phase-difference images; 

remove the non-zero background phase from the phase- 
difference images; and 

determine phase associated with flowing spins and phase 
associated with stationary spins. 

[Claim 2] 2. The MR scanner of claim 1 wherein the computer is 
further programmed to reconstruct a phase contrast image with contrast 
between flowing spins and stationary spins. 

[Claims] 3. The MR scanner of claim 1 wherein the computer is 
further programmed to: 

determine a degree of pulsatility of an object; 
compare the degree of pulsatility to a pulsatility threshold; and 
determine inclusion of the object in a background image from the 
comparison. 

[Claim 4] 4. The MR scanner of claim 3 wherein the computer is 
further programmed to exclude the object in the background image if the 
degree of pulsatility of the object exceeds the pulsatility threshold. 
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[Claim 5] 5. The MR scanner of claim 1 wherein the computer is further 
programmed to generate a mask image from those objects not pulsating at a 
temporal frequency above a frequency threshold. 

[Claim 6] 6. The MR scanner of claim 5 wherein the frequency 
threshold is 0 Hertz. 

[Claim 7] 7. The MR scanner of claim 6 wherein the computer is 
further programmed to generate a temporal frequency power spectrum on a 
per pixel basis, and if a given fraction of the temporal frequency power 
spectrum for a given pixel exceeds the frequency threshold, exclude the given 
pixel from inclusion in the mask image. 

[Claims] 8. The MR scanner of claim 7 wherein the computer is 
further programmed to exclude a pixel from inclusion in the mask image if 
more than 1 5 percent of pixel energy is at or above the frequency threshold. 

[Claim 9] 9. The MR scanner of claim 5 wherein the computer is 
further programmed to subtract the mask image from a phase-difference 
image to remove that portion of the phase-difference image corresponding to 
background phase. 

[Claim 10] 10. The MR scanner of claim 1 wherein the gradient coils are 
further configured to apply flow-encoding gradients with a non-zero first 
moment to encode spins with a phase proportional to their velocity. 



Page 1 9 of 30 



[Claim 11] 11. The MR scanner of claim 9 wherein the computer is 
further programmed to subtract background phase from the phase-difference 
image for each phase of a cardiac cycle. 

[Claim 12] 12. The IVIR scanner of claim 9 wherein the computer is 
further programmed to subtract background phase from the phase-difference 
image along a plurality of flow directions. 

[Claim 13] 13. The IVIR scanner of claim 3 wherein the computer is 
further programmed to determine inclusion of the object in a background 
image on a per slice basis. 

[Claim 14] 1 4. A method of phase correction in flow analysis MR imaging 
comprising the steps of: 

determining temporal frequency components for a given pixel of a 
phase-difference image; 

generating a power spectrum of the given pixel from the temporal 
frequency components; 

determining percentage of the power spectrum within a given energy 
range; and 

including the pixel in a mask image if a given percentage of the power 
spectrum for the pixel is at a given energy. 

[Claim 1 5] 15. The method of claim 1 4 further comprising the step of 
including the pixel in the mask image if less than 1 5 percent of the power 
spectrum for the pixel is at the given energy. 

[Claim 16] 16. The method of claim 14 wherein the given energy is DC. 
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[Claim 1 7] 17. The method of claim 1 4 further comprising the step of 
applying the masl< image to the phase-difference image and generate a set of 
masked phase-difference pixels. 

[Claim 18] 18. The method of claim 1 7 further comprising the step of 
determining a set of fit coefficients of a representative polynomial expression 
from the set of masked phase-difference pixels. 

[Claim 19] 19. The method of claim 1 8 further comprising the step of 
determining background phase in the phase-difference image from the set of 
fit coefficients. 

[Claim 20] 20. The method of claim 1 9 further comprising the step of 
subtracting background phase from the phrase-difference image and 
generating a corrected phase-difference image. 

[Claim 21] 21. The method of claim 20 further comprising the step of 
generating a corrected phase-difference image for each slice of a plurality of 
imaged slices. 

[Claim 22] 22. The method of claim 20 further comprising the step of 
subtracting background phase from the phase-difference image on at least 
one of a per flow direction and a per cardiac phase basis. 

[Claim 23] 23. The method of claim 1 8 wherein the representative 
polynomial expression is a first order expression. 

[Claim 24] 24. The method of claim 1 4 further comprising the steps 
summing the temporal frequency component for a first harmonic and the 
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temporal frequency component for a second harmonic of the power spectrum 
for the given pixel and normalizing the sum to a maximum pixel value of the 
phase-difference image. 

[Claim 2 5] 25. The method of claim 24 further comprising the steps of 
comparing the value of normalization for the given pixel to a threshold and 
determining inclusion of the given pixel in the mask image from the 
comparison. 

[Claim 26] 26. The method of claim 25 further comprising the step of 
applying a pre-selected threshold to determine inclusion of the given pixel in 
the mask image. 

[Claim 27] 27. The method of claim 25 wherein the pre-selected 
threshold is set to 0.04 such that those pixels having a value of normalization 
above 0.04 are excluded from inclusion in the mask image. 

[Claim 28] 28. A computer readable storage medium having a computer 
program stored thereon and representing a set of instructions that when 
executed by a computer causes the computer to: 

(A) determine a pulsatility factor for each pixel of a slice of a phase- 
difference image; 

(B) generate a mask image of the slice from only those pixels having a 
given pulsatility factor; 

(C) apply the mask image to the phase-difference image for the slice to 
mask the phase-difference image along the slice; 

(D) repeat (A)-(C) for each slice of the phase-difference image; and 
(E) output a corrected phase-difference image with substantial 
subtraction of background phase. 
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[Claim 29] 29. The computer readable storage medium of claim 28 
wherein the set of instructions further causes the computer to apply the mask 
for each cardiac phase of a measurement period. 

[Claim 30] 30. The computer readable storage medium of claim 28 
wherein the set of instructions further causes the computer to apply the mask 
along each direction of flow in the slice. 

[Claim 31] 31. The computer readable storage medium of claim 28 
wherein the set of instructions further causes the computer to determine 
background phase in the phase-difference image with perturbation as a 
function of inflowing spins along the slice. 

[Claim 32] 32. The computer readable storage medium of claim 31 
wherein the set of instructions further causes the computer to calculate 
quantitative flow measurements to access flow velocities in vivo. 

[Claim 33] 33. The computer readable storage medium of claim 28 
wherein the set of instructions further causes the computer to low pass filter 
the phase-difference image and Fourier transform the filtered phase- 
difference image to determine temporal frequency components for each pixel. 

[Claim 34] 34. The computer readable storage medium of claim 33 
wherein the set of instructions further causes the computer to sum first and 
second harmonics of a power spectrum for each pixel, normalize the sum to 
maximum pixel value in the phase-difference image, and compare the 
normalization value of each pixel to a threshold. 

[Claim 35] 35. The computer readable storage medium of claim 34 
wherein the set of instructions further causes the computer to use a pixel in 
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generation of the mask image if the pixel has a normalization value less than 
the threshold. 

[Claim 36] 36. The computer readable storage medium of claim 28 
wherein the phase contrast image is a speed image of encoded flow along a 
plurality of flow directions and wherein the set of instructions further causes 
the computer to determine velocity components along the plurality of flow 
directions. 

[Claim 37] 37. The computer readable storage medium of claim 36 
wherein the set of instructions further causes the computer to generate an 
individual power spectrum for the velocity components and determine 
inclusion/exclusion of a pixel in the mask image from the individual power 
spectrum for the pixel. 

[Claim 38] 38. The computer readable storage medium of claim 37 
wherein the set of instructions further causes the computer to include the 
pixel in the mask image if a given percentage of the power spectrum for the 
pixel is at 0 Hertz. 
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